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Abstract We consider the harmonic chain of oscillators with self-consistent 
stochastic reservoirs and give a new proof for the finitude of its thermal 
conductivity in the steady state. The approach, with involves an integral 
representation for the correlations (heat flow) and a perturbative analysis, is 
quite general and extendable to the study of anharmonic systems. 
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In nonequilibrium statistical physics, the analytical derivation of the macro- 
scopic laws of thermodynamics for microscopic models of interacting particles 
is a challenge to theorists. In particular, even in ID context, it is still unknown 
the rigorous derivation of Fourier's law of heat conduction which states that 
the heat flow is proportional to the gradient of the temperature [T]. Many 
works, almost all of them by means of computer simulations [2] and with 
conflicting results have been devoted to the theme since the pioneering rig- 
orous study on the harmonic chain of oscillators with thermal baths at the 
boundaries [3J, a model that does not obey the Fourier's law. Recently, in 
the scenario of analytical studies, the harmonic chain of oscillators has been 
revisited, but in the case of each site coupled to a stochastic reservoir: it is 
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proved that the Fourier's law holds in such case [4]; of course, if we turn off 
the coupling with the inner reservoirs, the heat conductivity diverges and 
the Fourier's law does not hold anymore, as previously shown in [3]. As the 
procedure presented in [3] was completely dependent on the linearity of the 
dynamics (i.e., on the harmonicity of the potential), some of the present au- 
thors have proposed a different approach to study the chain of oscillators 
with a reservoir at each site in order to analyze the problem of anharmonic 
potentials and infer their behavior with thermal baths at the boundaries 
only J5][6l[7]. The proposed approach involves an integral representation for 
the correlations (related to the heat flow), whose analysis is carried out by 
means of a perturbative computation. 

In the present letter, as a rigorous support for our previous theoretical 
works (which involve, as said, anharmonic interactions), we turn to the sim- 
pler case of the harmonic chain with a reservoir at each site and give a new 
proof for the finitude of the thermal conductivity (that shall be extendable 
to the anharmonic case), using our approach and a perturbative analysis: 
the expression for the conductivity presented here coincide with that one 
described in [4|. 

Let us introduce the model. We consider a lattice system with unbounded 
scalar variables in a space A C Z and with stochastic heat bath at each site. 
Precisely, we take a system of N oscillators with the Hamiltonian 

N N 

H ^ 1) = E 2 $ + M ® + qjJjm > (1) 

where Jy = Jji and Jij = fQl — In the study of harmonic chains it is 
usually taken j and I nearest neighbors and Jj,j+i = J for any j. In this 
case, the interparticle interations may be also written as (assuming Dirichlet 
boundary conditions qo = = qN+i) 



N 

"£V( qj -q j+i ) , V = 




3=1 

after adjustments in M, the coefficient of q|. We consider the Langevin time 
evolution given by the stochastic differential equations 

dqj = pjdt, 

dpj = —£-dt - frjdt + i^dBj, j = 1, 2, • • -N, (3) 
aqj 

where Bj are independent Wiener processes (i.e., rjj = dBj/dt are indepen- 
dent white noises) £ is the coupling between the site j and its heat bath; 
jj = 2£7j, where Tj is the temperature of the j-th reservoir. As usual, we 
define the energy of a single oscillator as 

Hj(q,p) - \v) + qj + Ufa) + \^ Vfa - «), (4) 
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where the expressions for U and V follow from eq. ([T]) and from J^j Hj = H . 
Then for the energy current we have 

where (•) means the expectation with respect to the noise distribution; Rj 
denotes the energy flux between the j-th reservoir and the j-th site 

(R j )=Z{T j -(p 2 j )); (6) 

and the energy current in the chain is 

l>3 



Pj+Pi 
2 ' 

Pi + Pj 



(7) 



As well known, the stationary state is characterized by (dHi/dt) — 0; and, for 
physical reasons, it is interesting to consider the "self-consistent condition" 
given by (Rj) = in the steady state. For the linear dynamics, the existence 
and convergence to the stationary state as t — > oo are old solved problems, 
see e.g. [8]. 

Lets us state our main theorem. 



Theorem 1 For the harmonic chain of oscillators with reservoirs at each 
site (1-3), in the case of nearest neighbor interactions, i.e. Jji = J((5; J+ i + 
5i t j—i), with J < Jo for some small Jo, the Fourier's law holds 

T = lim (TjJ) = -^^{Tn - Ti), (8) 

t-*oo JS — 1 

with the "self consistent" condition lim^oo (Rj) — for the inner sites j, 
and with the heat conductivity x given by 

Remark 1 A theorem establishing the finitude of the thermal conductivity for 
the considered model has been already presented in [3], as already said, and 
with the nearest neighbor interparticle potential not necessarily small. Our 
aim in this letter, as emphasized, is not to give a second proof but to present 
a more general approach which shall extend to the anharmonic chains. 

Remark 2 Following the steps to be presented ahead, we may also study the 
heat flow for weak interparticle interactions beyond nearest neighbor: e.g. for 
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Now we describe our approach. It is usefull to introduce the phase-space 
vector (f> = (q,p) and write the dynamics as 

4> = -A(f> + at], (10) 
where A = A Q + J and a are 2N x 2N matrices 

^=UrV=(;s)-=(«v^)< <"> 

where / is the unit N x N matrix, J is the N x N matrix for the interparticle 
interactions, and A4, r and T are the diagonal N x N matrices: Mji = M5ji, 
Fji = £6ji, T = TjSji. And rj are independent white noises. In what follows we 
will use the index notation: i for index values in the set [iV+1, 7V+2, • • • , 2N]; 
j for values in the set [1, 2, • • • , N] and k for values in [1, 2, • • • , 27V]. We will 
also omit obvious sums over repeated indices. 

Our strategy is first to consider (fT0|) above with J = 0, i.e., a system 
with isolated sites (without interactions among them). In a second step we 
introduce the interparticle interaction by using the Girsanov theorem, and 
then we calculate the heat flow in a perturbative computation. 

We have 

Lemma 1 The solution of \ 10\) with J = is the Ornstein- Uhlenbeck Gaus- 
sian process 

(P{t) = e- tA 'V(0) + / dse- {t - s)AO an{s), (12) 

where, for the simple case of 0(0) = 0, the covariance of the process evolves 
as 

utAM * -ru \ S e- (t - s)A V( s , s ), t>s, 
( ^( s))o = C (M) = | c(M)e _ (s _ t)AO T ; ^ 

C(t,t)= f dse- sA °a 2 e- sA ° T 
Jo 

Proof It is a simple exercise of stochastic differential equations (see e.g. [9], 
p. 74 )□. 

In the simple case of J = 0, as t — > oo we have a convergence to the Gaussian 
distribution with covariance 

C = J™ dse- sA °a 2 e- sA ° T = ( f ° ) , (15) 

where as said , T is a diagonal N x N matrix with %j — TiSij (see e.g. [5]). 

The solution of (fT0|) with the interparticle potential will be derived from 
the particular case of J = by using the Girsanov theorem which establishes 
a measure p for the complete process in terms of the measure ]ic obtained 
for ^7 = 0. Precisely, for the two-point function we have 



(13) 
(14) 
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Lemma 2 The two-point functions for the complete process 110\) can be writ- 
ten as 



((Pu(ti)<p m (h)) = / 4>u(ti)4> m (t2)Z(t)d{ic/norm., ti,t 2 < t, 



(16) 



where (f> is the solution (given by Lemma 1) of the process with J = 0, and tp 
is the solution for the complete process C is given by Iil3\14\ ), and the 

corrective factor is 



Z(£)=cxp^y u-dB-^J u 2 ds 



V 2 T J 

7» u i = -Ji-Nj<pj- 



(17) 



Proof For the harmonic potential, the process is an Ito diffusion, and so, the 
proof is also direct: see e.g. theorem 8.6.8 in [S]. □ 

Remark 3 In the case of the nonlinear (anharmonic) dynamical system prob- 
lem we may introduce the anharmonic interactions (on-site and/or interpar- 
ticle potentials) still by using the Girsanov theorem: e.g., for bounded po- 
tentials or initial processes with continuous realization, the use of Novikov 
condition makes the procedure straightforward; see an example of other ma- 
nipulations in a similar use of the Girsanov theorem [lOj . 

Turning to the Z(t) expression above, we have 
mdUi = 7 4 mjt dBi 

= 1 ; 1/2 u l (#i + A° k ct> k dt) 
= - 1 ; 1,2 j l ^ 3 {d4> l + A° lk 4> k dt). 

Using the Ito formula we get 

';/,,<.v/o, = -dF - f^faJijA^fadt, 
F{4>)=^ ( j )i Ji j <f>o- 



And so, 



Z{t) = exp / u-dB 



u ds 



exp [-F(0(t)) + F(0(O))] exp 



W(<p(s))ds 



W(4>(s)) = l^U^JiA^ + ^Af^^J^ + 



(18) 
(19) 

(20) 
(21) 



+ ^j'( s )Jj'i7i JijM 3 )- 

To analyze the heat flow in the steady state (related to lim^oo (ip u (t)ip v (t))), 
we note that we may write 



exp (-tA°) = e-*l cosh(ip) 1(q°j 



tanh(fp) 
P 



2 

-M 



I 



(22) 
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/ 2 \ 1 / 2 

p = ( f|J — Mj .It may be shown by e.g., diagonalizing A (or see the 
appendix of [3]). We also note that 



C(t, s) = exp (-{t - s)A°) C + G(exp (-(t + s)£/2)), 



(23) 



where the effects in the correlation function formula of the second term on 
the right-hand side of the equation above vanishes in the limit t — > oo; C is 
the covariance (fT5|) . 

Now, to compute the two-point correlation function (|16p (and so, the 
heat flow), we expand the exponential which gives Z(t) above (|2QII2ip in a 
power series: expLY] = '}2^ =Q X n /n\ 1 and calculate the connected Feynman 
graphs for (|16[) fdue to the normalization factor we stay with the connected 
graphs only) using the Wick theorem. As we have quadratic terms in (f> in 
the exponent (see F and W above), there is no countable problem with the 
series; roughly, 



b 2n 



—rd/ic 



(2n-l)M 



3 E <££: 



■C) 



fe=l 



■C\=2 n \C---C\. (24) 



In short, to control the expansion we only need a bound for the n convolutions 
of the covariance C such as c" J™, where c is some constant. Thus at least for 
small J, this simple analysis will give us the convergence of the perturbative 
series. To get such bound, we use the formulas (|22ti23[) for C, the lemma below 
and that 



||C|| < - e- 2 "*), 
\e- tA °\\ <c 2 e- Qt , 
P°ll<c 3 , 



where a = min{£/2, M/£} (|| • || means a matrix bound on M2NX2N 
Lemma 3 Let I t be 



It = 



f e -«|t-sil e -«| S i-s2| ...e'^-^dsf-dsn, a>0, (25) 
Jo 



then, limt—joo I t < (c a ) n , where c a does not depend on n. 
Proof For f{x) = e - "' 21 ', we have 
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And 



DC- 



f(t - «i)/(si - s 2 ) ■ ■ ■ f(s n - t)dsi - --dsn 



= -/*/* •••*/(<)) 
- v „ ' 



n+1 



1 Z" 00 
= -(2tt)("- 1 )/ 2 / (/(?))"+>, (27) 

where * means the convolution, and we have used Parseval's theorem in the 
last equality above. Hence, 

lim I, < - ( 2.)<-^ y (-^) * < (c). (28) 



Thus, using the formulas (|22j|23p for the covariance and the lemma [3] above, 
we may bound the terms order J 2 and up in the two point function (|16p by 

oo 

5Z(c'J) n < c"J 2 , (29) 

n=2 

where c' and c" are some constants (i.e., they do not depend on n). 

To obtain the result of theorem [T] we still need to calculate, in detail, the 
two-point correlation function up to first order in J. Carrying out the com- 
putation (simple Gaussian integrations), for (ip u fv) = lim^oo { L Pu{t) L Pv(t)) , 
we obtain 



ue [N + !,■■■ ,2N], 
v £ [1, • • • , N], (30) 
T U _ N 5 U , V , for u, v e [N + 1 • • • , 2N] . 



hp \ — J 2^M \^v+N,u-nT u -N — Ju.vTv] , for 



Hence, from ([7]) 



F j ^=2_jJ j+N ,r{<Pj-<Pr)— ^ , re[l,---iV], (31) 



r>j 

and we have (for (tp u if v ) up to first order in J) 



<^>=E%fr^- T ;)' 02) 



or, for nearest neighbor interactions, 

J; 

2£_M 



^i-i+i = <^> = (Jj ™ +l)2 (Tj+i - T 3 ) . (33) 
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The first order perturbation computation (j30|) above still gives us 

lim (R 3 (t))=0. 

t—>oo 

Thus, the steady state condition (dHi/dt) — leads to 

•^1^2 = ^2^3 = ■ ■ ■ = Tn-i^n = J~ ■ (34) 

And so, for the simpler case of the same interactions between any two nearest 
neighbors sites, i.e., J = Jj+N,j+i for any j, it follows that (for (ip u <Pv) up 
to first order in J), we have 

Tn — Ti J 2 

^ =x ^r' x =m- (35) 

From ([31"]) and (|29|) , it is easy to see that considering the remaining terms 
(higher order in J), we get as J x cJ 2 = cJ 3 in \ (the temperature terms 
Tj + i — Tj may be extract from the perturbative analysis following the prod- 
ucts of C, the covariance (fl"5]) in C which appears in the Gaussian interaction, 
and the terms 7$ coming from the expansion of Z(t)). In short, theorem [T] 
holds. 

Remark 4 In the case of anharmonic terms in W and F (|19H21| . a naive 
expansion fails, but we expect to develop a perturbative approach by using 
some cluster expansion as well known in field theory and statistical mechan- 
ics. For the simpler case of nonconservative nonlinear stochastic dynamical 
model (describing a system in contact with thermal reservoirs at the same 
temperature, and so going to equilibrium), a convergent cluster expansion is 
presented in [10] , and the decay of the four-point function is investigated in 
[12] and [11] : the complete and rigorous result [11] adds only small corrections 
to the first order perturbative calculation [T2] , 
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